o 



o 
o 



s 



Wall emission in circumbinary disks: 
the case of CoKu Tau/4 

Erick Nagel^'^, Paola D'Alessio^, Nuria Calvet^, Catherine Espaillat^'^, Ben Sargent^, Jesiis 

Hernandez^, and William J. Forrest*^ 

Centra de Radioastronomia y Astrofisica, UNAM, Morelia, Michoacdn, Mexico 58089 

Departamento de Astronomia, Universidad de Guanajuato, Guanajuato, Gto, Mexico 36240 

Department of Astronomy, University of Michigan, Ann Arbor, MI 48109 



Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-78, Cambridge, MA, 
PLh' 02138 

pq; 

^ ! NSF Astronomy & Astrophysics Postdoctoral Fellow 

Q . Space Telescope Science Institute, Baltimore, MD 21218 

^ ■ Centro de Investigaciones de Astronomia, Merida 5101- A, Venezuela 



Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 

>\ 
CN , erick@astro.ugto.mx 

m 
(N 

O 



ABSTRACT 



A few years ago, the mid-IR spectrum of a Weak Line T Tauri Star, CoKu 
Tau/4, was explained as emission from the inner wall of a circumstellar disk, with 
k> ! the inner disk truncated at ~ 10 AU. Based on the SED shape and the assumption 

that it was produced by a single star and its disk, CoKu Tau/4 was classified as 
a prototypical transitional disk, with a clean inner hole possibly carved out by a 
planet, some other orbiting body, or by photodissociation. However, recently it 
has been discovered that CoKu Tau/4 is a close binary system. This implies that 
the observed mid-IR SED is probably produced by the circumbinary disk. The 
aim of the present paper is to model the SED of CoKu Tau/4 as arising from the 
inner wall of a circumbinary disk, with parameters constrained by what is known 
about the central stars and by a dynamical model for the interaction between 
these stars and their surrounding disk. We lack a physical prescription for the 
shape of the wall, thus, here we use a simplified and unrealistic assumption: the 
wall is vertical. In order to fit the Spitzer IRS SED, the binary orbit should 
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be almost circular, implying a small mid-IR variability (10 %) related to the 
variable distances of the stars to the inner wall of the circumbinary disk. In the 
context of the present model, higher eccentricities would imply that the stars 
are farther from the wall, the latter being too cold to explain the observed SED. 
Our models suggest that the inner wall of CoKu Tau/4 is located at 1.7a, where 
a is the semi-major axis of the binary system (a ~ 8AU). A small amount of 
optically thin dust in the hole (< 0.01 lunar masses) helps to improve the fit to 
the 10 yum silicate band. Also, we find that water ice should be absent or have 
a very small abundance (a dust to gas mass ratio < 5.6 x 10^^). In general, for 
a binary system with eccentricity e > 0, the model predicts mid-IR variability 
with periods similar to orbital timescales, assuming that thermal equilibrium is 
reached instantaneously. 

Subject headings: circumstellar matter - infrared: stars - stars:pre-main-sequence 



1. Introduction 



Until recently, the Weak Line T Tauri star, Co Ku Tau/4 (|Cohen fc Kuhilll979l ). was 



known as a remarkable exam ple of a transitional dis k (JForrest et al. 



2004). Its mid-IR SED 



(IRS Spitzer) was modeled by lD'Alessio et al.l ( 120051 ). hereafter D05, taking an isolated star, 
and they concluded that the disk surrounding CoKu Tau/4 should have a clean inner hole, 
with a radius about 10 AU. In this model, the emergent flux in the mid-IR spectral range 
was produced by the inner wall of the truncated disk. 

Previously, it had bee n thought that the inner hole of CoKu Tau/4 was due to either 



a recently formed planet (Quillen et al.l |2004J ) or by the UV switch model ( IClarke et al. 



200 ll : [Alexander fc Armitaed 120071) due to i ts un detectable accretion signatures. However, 
observations described by [Ireland fc Kraud (120081 ) show conclusively that CoKu Tau/4 is a 
binary system, implying that its inner hole is probably produced by dynamical interactions 
between the binary system and the surrounding gas. 

For a binary system fo rmed by stars of sirn i lar rn a sses and using the pre dictions of 
Artymowicz &: Lubowl (Il994| ) and lBeust fc Dutreyl (|2005[ ) . Ilreland &: Kraua ( 120081 ) estimated 
that the size of the hole in the circumbinary disk (CB hereafter) should be between 13 and 
16 AU, for different eccentricities of the binary orbits. The estimated hole is larger than the 
one found in the D05 transitional disk model, and being further away, one would naively 
expect it to be cooler. However, the stars are not at the geometrical center of the disk: 
each point in the visible part of the wall is at a different distance from the stars. Therefore, 
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the net effect of botli stars in tlie fieating of tlie wall should be calculated for a particular 
configuration before comparison with the observed SED. 

This paper has several goals. The first is to find a simplified structure of the inner 
edge of the CB disk around the CoKU Tau/4 binary system that consistently reproduces 
the observed spectral energy distribution (SED). The second is to determine the properties 
of the dust particles of that edge. Spitzer ob servations probe mos tly the upper layers of 



the inner 10-20 AU regions of typical disks (JD'Alessio et al.M2006l ). while at the present 
resolution of submillimeter and millimeter interferometers, the main contribution comes 
from the midplane at larger distances. In contrast, circumbinary and transitional disks are 
truncated at a certain radius, making them a natural laboratory to study the dust near 
the midplane in the inner disk region, inaccessible by other means. In § [2] we describe the 
physics of the models, and in § [3] we present the fiducial model adopted as a reference for 
later comparisons; in § |4] we study the effects of changing parameters in the SED of the 
fiducial model; in § [5] we apply our methods to the case of CoKu Tau/4, and in § [6] we 
discuss our results and give conclusions. 



2. Description of the model 

Our basic model consists of a CB disk, from which only its inner wall contributes to 
the SED in the mid-IR spectral range (5 — 40/im), surrounding two stars that contribute 
mainly to the near-IR flux. The main difference between the present model and the one 
described in DOS is that here the inner wall is irradiated by two stars and each of them has 
a different distance to any surface element of the wall since neither of the stars is located at 
the geometrical center of the CB disk. Moreover, since the stars are moving in their orbits, 
their distance to each portion of the inner wall changes with time. Assuming that the wall 
reaches thermal equilibrium in a timescale shorter than the orbital timescale of the binary 
system (see the Appendix for the evaluation of this assumption), the temperature of each 
portion of the wall and the resulting SED are periodically variable with its orbital timescale. 

The fact that we are describing a CB disk implies that a whole set of new parameters 
should be included into the model, in comparison to the model proposed by DOS. At first 
sight, it might look like this increasing degree of freedom makes the problem degenerate. 
However, new restrictions allow one to break part of this degeneracy, as we will see in the 
following sections. 

Figure [1] shows the geometrical parameters that describe the spatial configuration of the 
system, namely, the semi-major axis a, eccentricity e and the angle 0, that locates one of 
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the stars in an elliptical orbit with the other one at the focus. Figure [T] shows the two angles 
that are required to characterize this configuration in 3D space: i, the inclination between 
the disk axis and the line of sight (or, between the disk plane and the plane of the sky) and 
a, which gives the axis in the disk plane around which the disk is inclined. This line passes 
through the center of the elliptical orbit and the angle a is measured between this line and 
the major axis line. Note that first the system is described in the orbital plane with a, e and 
(f) (see Figured]). Then, the major axis line is identified, from which an angle a gives the line 
around which the system is inclined an angle i. An inclined system is shown in Figure HI 
For these parameters, the origin of the reference system is in the center of the ellipse. Notice 
that changing a not only changes the orientation of the visible part of the wall, but also 
modifies the distances from the stars to that part of the wall (§ 12. ip . On the other hand, a 
variation in i changes the visible area of the wall (§ 12. ip . In the following, the angles and 
a are given in units of vr radians. 

The inner wall or inner boundary of the CB disk is characterized by its radius Rcb, and 
its vertical height, h. Notice that in this system of reference, while the second star describes 
the elliptical orbit, the center of mass of the binary system and also the center of the hole in 
the CB disk move (both po ints are the same for e = but not for other values of eccentricity, 
e.g., iPichardo et al.l (|2005[ ) and lPichardo et al.l ( 120081 )). Thus, during a binary system orbit. 



the distance between the two stars and the wall changes. For the following, the origin of our 
system of reference is the center of the CB disk. 

Also, it is important to specify the dust composition of the material in the hot atmo- 
sphere of the wall. We assume that, in general, the dust is a mixture of silicates, troilite, 
graphite, amorphous car bon and water i ce. T he silicate (olivine or pyroxene) and graphite 
abundances are given by iDraine fc Led (Il984[). The maxi mum abundance of water ice and 



the abundance of troilite is taken from lPoUack et al.l (J1994J ). We change the dust composition 
and abundances in §|l]and adjust them by comparing the synthetic and the observed SEDs 
in §0 

A constraint inherent to a two-star model is given by the dynamical interaction of the 
binary system with a disk. This interaction clears an inner hole, thus, in this way the 
radius Rcb is determined by the stars' masses and the eccentricity and semi-major axis of 
their orbits. In contrast, in the case of one star (DOS) Rwaii is a free parameter. In order 
to reproduce the observed SED, Rwaii was adjusted to obtain the appropriate temperature 
structure. 



The size o f the hole (Rrh) can b e calcu l ated numericallv (IRudak fc Paczvnskil ( llQBll ). 
Pichardo et al.l ( 120051 ) JPichardo et al.l (20081 ). JArtvmowicz fc Lubowl (19941)') o r ana lytically 
for a circular binary (JNagel et al.l l2008l ). However, lArtymowicz fc Lubowl (Il994[ ) take a 
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hydro dynamical approach with the inclusion of viscosity, in addition to the gravitational 
interaction of the stars and the disk material. The presence of viscosity allows the material 
to move to the region where there are no stable orbits around both stars. If a mass par- 
ticle passes a certain point, then the only possibility for it is to move to one of the stars. 
In this way, streams of material originating in the CB disk connect to the circumstellar 
disks (these are disks inside the hole surrounding ea ch star), acting as a feeding mec h anisrn 
( JArtymowicz fc Lubowlll996l : iGiinther fc Kleyll2002l ). Observations by IJensen et al.l (|2007| ) 
suggest the presence of st reams from the CB di sk towards the binary system UZ Tau E. 
Also, for the binary St 34, iHartmann et al.l ( 120051 ) argue that the presence of dust in the dy- 
namically cleared hol e is consistent with having some accretion flow onto the central binary. 
Giinther et al.l ( 120041 ). described the emission of DQ Tau and AK Sco using a hydrodynami- 
cal simulation, where all the disks and the streams connecting them are tak en into account. 
Also, binaries of very low mass are expected to be fed by a circumbinary disk (JRobberto et al. 
20081 ). From all these studies, we can argue that the truncatio n radius Rrh strongly d e pends 
on the viscosity in a non trivial way. For an e > system, lArtymowicz fc Lubowl (119941 ) 
show clear differences between Rcb for changes in the Reynolds number, which depends on 
the vis cosity. However, in order to take a more realistic model, we decided to use the re- 
sults of lArtymowicz &: Lubowl (1l994l ) taken from Figure 3, arguing that Rcb does not depend 
strongly on the stars mass ratio, and a Reynolds number of 10^. 

An important ingredient for the calculation of the radiation flux from the stars impinging 
on the circumbinary wall is the distance between each star and a specific point on the wall, 
Ri,2 (see Figures H] and |2]) . In order to calculate these distances, the first step is to give the 
locations of the stars in the configuration where one star is at the focus and the other is on 
a point on its corresponding ellipse. The second step is to place a circular CB disk with its 
origin at th e point corresponding to the eccentricity of the binary system (see lPichardo et al. 



2005 



2008r). e.g., the c enter of mass of the binary system if e = 0, with a radius given by 



Artymowicz fc Lubowl (Il994j ). Finally, the third step is to calculate the distances between 
the coordinates of the stars and each surface element of the wall given in this reference 
system (see Figures [1] and |2]) . 



2.1. Formulation of the wall emission problem 



The inner wall of the CB disk is assumed to be vertical, i.e., the vector normal to its 
surface is parallel to the disk midplane. In order to relax the assumption of verticallity one 
requires a way to describe the shape of the wall, which is not possible to the best of our 
knowledge, because the observations and simulations do not have the necessary resolution. 
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More complicated shapes have been assumed for the dust to gas transition at the dust 
destruction radius in th e inner disk by taking into account the dependenc e of the subhmation 



temp erature or density (jlsella fc Nattall2005l ) or the maximum grain size (JTannirkulam et al. 



20071 ). However, the edge of the CB disk is a physical wall created by tidal interactions 
with the binary system, and the vertical density profile may be very different from that in 
standard disks. Since no detailed models exist to guide us, we have adopted the simplest 
geometry: a vertical wall. The wall is also assumed to be circular for simplici ty, however, non- 
axisym metr ic shapes are expected b ased o n hydrodynamical simulations of iGiinther fc Kiev 
(I2OO2I ) and lArtymowicz fc Lubowl (119961 ). In Section 14.11 we show that the SED for an 
eccentric CB disk does not change appreciably with respect to the circular case, thus the 
CoKu Tau/4 modeling is restricted to the latter. 

Only a fraction of the wall's surface can be detected by the observer, depending on the 
wall radius, height and the inclination angle i. There are other portions of the wall that 
are occulted from the observer by absorption due to the disk itself. The "visible" surface, 
projected in the plane of the sky (see Figure [2]), is divided in small surface elements (pixels). 
Each one is located at a given height relative to the disk midplane, Y and at a given distance, 
X, from the disk center. In order to describe each pixel position, we use cartesian coordinates 
(see DOS). 

Each point in the wall has a different temperature distribution, from the surface toward 
larger radii (deeper inside the dis k) . In order t o estim ate t his temperature dist ribution we 
use a similar treatment given by ICalvet et al.l (1199 if ) and ICalvet et al.l (119921 ). However, 
in the present case each pixel is being irradiated by two stars, neither of which is at the 
center of the system. We use the fact that the luminosity from the wall is much less than 
the luminosity from the stars in order to neglect its contribution to the heating of the wall. 
Thus, at each pixel we consider the superposition of two radiation fields, characterized by two 
different fluxes and mean intensities, and, in principle, by two distinct incidence directions. 
The incident flux from each star is F^j = {L^j/A7iR'j)fi with each star at a different distance 
from the wall, given by Rj, /x is the cosine of the incidence angle with respect to the normal 
to the wall surface. Again, for the sake of simplicity, we assume that each pixel is a plane 
parallel atmosphere and the radiation from each star arrives in a single beam, with a normal 
incidence angle {fi = 1). This assumption is valid only if the minimum distance between 
the stars and the wall is much larger than the stellar radii. In the case of CoKu Tau/4, the 
minimum distance between one star and the wall is 1.2 a ~ 9AU, so the later assumption is 
justified for this particular object. 



Following ICalvet et al.l (Il99ll ) and DOS, the radiation field arriving at each pixel is 
then separated into three components: "stellar 1", "stellar 2" and "disk", each one with a 
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particular wavelength range where most of its emission peaks. Defining the optical depth at 
the "disk" wavelength range in the radial direction, measured from the wall surface towards 
larger radii, the net outward Eddington fiux at the disk frequency range, Hd-, is given by 

2 p 

Haira) = J2 ^"i[(^ + C'lje-"^-^'' + C2,,e-^^*-''], (1) 

i=i 

where aj is the mean true absorption coefficient, Ksj, divided by the mean total disk opac- 
ity, Xd (calculated including true absorptions and scattering), evaluated using the effective 
temperature of each star (given by j = 1 or 2), and (3j = {^ajY^"^. Here, qj{= Xs,j/Xd), 
is the ratio between the mean total opacity at the stellar wavelength range and the mean 
total opacity at the disk range (which is the opacity used to calculate r^). The constants in 
equation ([T]) are given by equations 2 and 3 in DOS, but in here they are evaluated using the 
mean opacities corresponding to each star. 

All the mean opacities are calculated using the Planck function evaluated at the local 
or at the stellar effective temperature as the weighting function. The reason is that this kind 
of average is easy to calculate and the mean contributions of scattering and true absorption 
can be added together to give the total opacity, in contrast to a Rosseland kind of average 



(JD'Alessio et al.l Il998l ). Using this kind of mean opacity we obtain the correct value of 



temperature at Th. — >■ 0, b ecause this is adequate for the optically thin layers of the wall 



atmosphere (JMihalaslll978l ). 



The Eddington flux leads to the mean intensity at the disk frequency range, using the 
flrst momentum of the transfer equation and the usual boundary condition Jd{0) = 2Hd{0). 

From the zeroth moment of the transfer equation, assuming isotropic scattering, the 
temperature can be written as, 

T{T,r = E "^-^(^i.. + ^2.e~^^- + C^,e-^^«^-), (2) 

with constants given in DOS. Another, more compact, way to write down equation (jj]) is 



T(r,)4 = Tl(r,)^ + T2(r,)^ (3) 

where Ti and T2 are the temperatures calculated as if only the star 1 or 2 are present, 
respectively, as a superposition solution. However, in spite of its apparent simplicity, this 



is really an implicit equation for the temperature since the mean opacities involved in the 
calculation of Ti and T2 depend on the temperature T. 



As in ICalvet et al.l ( 119911 ). ICalvet et al.l (Il992[ ) and DOS, we have assumed that the 
opacities are constant in the wall atmosphere, i.e., aj and Qj are independent of t^. If the 
opacity is dominated by dust, and the temperature distribution in the wall atmosphere is 
such that there is no sublimation of any dust ingredient at a particular optical depth, then 
it is valid to assume that aj is constant. Also, if the contrast in temperature between the 
uppermost and the deepest layers is not very high, Xd and Qj can be assumed to be constant 
too. Thus, in order to quantify the mean opacities, we use the temperature evaluated at 
Td = 0, i.e., the surface temperature T^aii- This is found solving the implicit equation ([2]), 
for a given wall radius Ret, a given configuration of the binary system, and the coordinates 
of each pixel in the wall, a pair {X, Y). 

The monochromatic emergent intensity from each pixel of the wall in the direction of 
the observer is calculated by integrating the transfer equation. We consider two components: 
the thermal emission and the scattering of stellar radiation. Finally, the SED is calculated 
by multiplying the emergent intensity by the solid angle subtended by each ("visible") pixel, 
as seen by the observer. 



2.2. Formulation of the hole emission problem 

As described in § [21 the size of the hole (Rcb) is given by the dynamical interaction 
between the disk and the stars. However, the full hydrodynamical treatment should include 
ingredients like the gas viscosity, that allows material from the CB disk to lose angular mo- 
mentum and fall into the hole. This results in a small quantity of material inside the CB disk 
hole, that can be treated as being optically thin. Avoiding the precise description of the dy- 
namical evolution of such particles, in this work we simplify the problem considering them dis- 
tributed in a region located between a radius R = Rmin and R = Rcb, uniformly distributed 
in the azimuthal angle in a cylindrical coordinate system, centered a t the center of mass of 
the b i nary system. Indeed, s imulations of the material inside the hole (jArtymowicz fc Lubow 
19961 : iGiinther Sz Kleyll2002l ) indicate that the material is non axisymmetric distributed and 
spans all radii less than Rcb- However, for a circular binary system with undetectable ac- 
cretion towards the stars, we can argue that most of the material is restricted to lie at 
radius larger tha n the outer Lagrangian point for the three body circular restricted problem 
(ISzebehelyl 119671 ) . The surface density is gi ven by a power law (S oc -R^ , p = 1) such that S 
increases closer to the wall, as is expected (jArtymowicz &: Lubowlll996l ). 
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The emergent flux from tlie dust in tlie liole depends on tlie dust temperature. Since we 
assume tlie hole is optically thin, then each dust grain is heated by absorption of geometrical 
diluted radiation from each star and cools down by emission of its own thermal radiation. 
This is expressed by the radiative equilibrium equation, written as 



K,Jldu= / K,B,{T{R,e))diy, (4) 

Jo 

where the mean intensity of the radiation field is approximated by 

j; ^ Wr{R^,R,r)B,{%^) + 1^2(i?2, i?*2)5.(T,2), (5) 

with the dilution factors. 



^ 2 



l-Wl-(^ 



,forj = l,2. (6) 



When the integrals are performed, the dust grain temperature is the solution of the 
implicit equation 

2 

T{R,ef = Y,W,Tt^[K,jKa{T{R,d))l (7) 

i=i 

where Kd{T{R, 0)) is the Planck mean opacity given by 

and K,s,j{T{R,6),T^j) is a Planck-like mean opacity, which can be written as 

.„(T(fl.»).r.,) = ii^!^(^l™^. (9) 

Here, we use the dust temperature to evaluate the monochromatic opacity and the 
stellar effective temperature (for each star, j = 1 or 2) to evaluate the Planck function in 
the average. The Planck-like mean opacity is the straightforward result from the solution of 
the radiative equilibrium equation. 
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The monochromatic opacities we use are calculated assuming a size distribution for the 
grains. In the present work, which is our first approximation to the problem, we do not 
calculate a particular temperature for each grain size, but a mean temperature for the whole 
size distribution. In this approximation, the emission is produced from small to large grains 
at the same temperature, however, for the full case, the emission from small and large grains 
will differ due to their different temperature. Thus, we assume that the differences in the 
emitted spectrum between the approximation and the real case are small. In § |4] we explore 
different ingredients for the dust in the hole. 



3. Fiducial model 



As described, this problem is too complex. Moreover, a full analysis of all the parame- 
ters involved is beyond our actual computational resources. Thus, we restrict the problem 
to variations around a fiducial model (El in table [2]), where some parameters are fixed. 
Observations of the stellar parameters and orbital parameters help restrict the set of free 
parameters. Comparison between the synthetic and observed SED then provides a range for 
the rest of the parameters. 

We adopt parameters relevant to CoKu Tau/4 for the fiducial model. First of all, the 
model requires the parameters that characterize the stars. What is ass umed for the stars 
is mostly based on the observations reported by [Ireland fc Kraud (120081 ) for CoKu Tau/4. 
The masses of the stars (M^i and M^2) have a ratio M^2/j^*i = 0.85 ± 0.05, and near-IR 
colors consistent with typical T Tauri spectral types. These stellar masses, combined with 
a, give the distance of each of the stars to the center of mass of the binary system. Also, 
assuming the stars are coeval, by specifying each of their masses, their ages and a stellar 
evolutionary model, one obtains the stellar radii (-R*i,2), their effective temperatures (T^i,2) 
and their luminosities (L^i^2)- 



In this paper we adopt the stellar evolutionary tracks from ISiess et al.l (|2000[ ). and the 
spectra are taken from synthetic stella r spectrum from libraries from A. Bruzual, which are 
described in lBruzual fc CharlotI ( 1l993l ). For practical purposes, we take the sum of the spec- 
tra of two appropriate stars for comparison with the obser ved spectra at la r ge fre quencies. 
This process results in two spectral types for both stars. [Ireland &: Kraud (120081 ) propose 
the stars have Ml. 5 ± 0.5 spectral types, co nsistent with an age of 4 Myrs. Here, we use the 
evolutionary tracks from ISiess et al.l ( 120001 ) given in O.IM© steps. Taking this restriction, 
we adopt MO and Ml spectral type stars, which are within the uncertainties (A. Kraus, 
personal communication). Table [1] shows the parameters for the stars. 
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Figure |3] shows the stellar spectra for the bin ary stars with an age of 4 Myrs, corrected 



using A^/ = 3 and four different extinction laws (IMoneti et a. 



20091: lMathi^ll990[ ). The spectra corrected with lPrainel J2003[ ) JMathisI Jl99o[ ). and 



2001 



Draind 12003 : 



McClure 



McClure 



(120091 ) are very close. Note that the extinction law in the near-IR does not introduce notice- 
able changes in the SED, but it affects the mid-IR spectral region (see Figure [3]). 

We note that from the spectrum there is no evidence of circumstellar disks, since the SED 
of CoKu Tau/4 at A < Sfim see ms to be completely described by photospheric emission of 
the stars (jlreland fc Krausll2008l ). Moreover, CoKu Tau/4 is classified a WTTS, with an Ha 
equivalent width consistent with negligible accretion from any possible circumstellar disks 
to the stars, with a rate constraine d to be M < 10~^^M(^/ year, which is an observational 
threshold for accretion, according to lMuzeroUe et al.l (120001 ) . Thus, in the fiducial model we 
do not consider the presence of circumstellar disks. 

We find that only a particular range in wall temperatures [Tu,aii{X, Y, r^)] produces the 
observed SED consistent with the observed SED of CoKu Tau/4. T^aii strongly depends on 
the range of distances of the stars to all the points in the wall. As noted above. Ret is no 
longer a free parameter and depends primarily on the eccentricity (e): if e increases, Rcb 
increases, too, and the i i both T^jaii and Fobs decrease. For the fiducial model, we adopt the 
Artymowicz fc Lubowl (119941 ) value for Rcb (= 1.7a) for a e = system for simplicity. 



In order to fix a and (p, we consider configurations with the highest emergent flux. 
For this to be done, we need to search between configurations along the whole orbit of the 
system. However, the angle is irrelevant because the orbit is circular, but the angles a 
and i have to be specified. The former defines what fraction of the wall is directed towards 
the observer and the latter, the projection of the wall in the plane of the sky. The stars are 
not illuminating the wall homogeneously and there are two intensity peaks in the positions 
(X,Y) closest to each of the stars. The only configuration in which both intensity maxima 
are observed is the case where the stars are located in the axis around which the system is 
inclined, = 0.5 and a = 0. Thus, we choose this configuration in an effort to reduce the 
number of free parameters and to get a reasonable model to start with. This also turns out to 
be the highest flux model. For the same reason, we take the minimum observed semi-major 
axis (a = 7A3AU) as the fiducial value. As in DOS, we choose cosi = 0.5 as a typical value. 
In the fiducial model the height of the wall is assumed to he h = 0.28a. 

Finally, for the fiducial model the dust composition consists of graphite, si licates, troilite 
and n o water ice, with dust to gas mass ratios (grap = 0.0025 and Ci( = 0.0034 (IDraine fc Lee 



19841 ). The troilite abund ance is Ctrni = 7.68 x 10~^ (IPollack et al.lll994l ). Graphite optical 
properties are taken from IWeingartner fc Draind (120011 ) for graphite, and the silicates are 
assumed to be pyroxenes {Mgo,s Feo,2 SiO^), with optical properties from iDorschner et al. 
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( 119951 ) ■ Also, we consider an interstellar medium dust size distribution (JMathis et al.lll977| ) 
with n{a) oc a"^'^, with minimum and maximum sizes amin = 0.005/im and ttmax = 0.25/im, 
respectively. 



4. Effect of parameters on SED 

To show the dependence of the SED on the input parameters we change one parameter 
at a time and compare the resultant SED to that of the fiducial model. We first consider 
the change of parameters in the wall, and then the change of parameters in the optically 
thin region inside the hole. We assume that the apparent separation of the stars on the 
plane of the sky is the s ame as the observed separation dap = 53.6 ± 0.5mas = 7.5 ± 0.07 AU 



(jlreland fc Krausll2008l ). This assumption implies a fixed semi- major axis, a. Note that the 
system configuration changes with time, thus, the dap will change, for a fixed a. 



4.1. Effects of wall parameters on SED 

In the modeling, the inner boundary of the CB disk is circular, which is the simplest 
assumption. However, in principle its shape can be non-axisymme tric. From the results 



of the hydrodynamical simulations of lArtymowicz fc Lubowl (119961 ) and iGiinther &: Kley 



(120021 ). we note that in the line that connects both stars, the disk material is approaching 
the outer lagrangian points. Thus, in order to quantify the difference between axisymmetric 
and non-axisymmetric models, we model the shape of the wall as an ellipse with its semiminor 
axis along this line. This is an easy way to get a non-axisymmetric wall, however, only a 
hydrodynamical simulation applied to systems like CoKu Tau/4 will allow us to accurately 
characterize the wall. Thus, the wall presented here is just one of many possibilities. 

Figure S] shows models with the parameters of the fiducial configuration, but with a 
semiminor axis for the CB disk of 6 = Rcb — (0.4, 0.3, 0.2, 0.1, 0.0)a. Note that the flux 
decreases for smaller b, which seems unphysical, due to the fact that the temperature of the 
wall closest to the stars is larger. However, the projected area perpendicular to the line of 
sight decreases in such a way that the second effect dominates the observed emission. 

The main result is that the changes are restricted to wavelengths longer than 12/im, 
thus, the silicates band is not modified. The changes are larger around 30/im, but this is 
only relevant for smaller b, which represents a highly eccentric disk. This example shows 
how an elliptical wall modifies the SED, however, only when observations and/or simulations 
solve the structure of the wall, we will be able to quantify this effect in the modeling of CoKu 
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Tau/4. 

Now, we consider changes of the binary system eccentricity, e, for = 0.5, a = 0, 
cosi = 0.5 and h = 0.28a. Figure shows the SED of the fiducial model and of models 



with e = 0., 0.2, 0.4, and 0.6 calculated with Rcb from lArtymowicz fc Lubowl (119941 ). As e 
increases, so does Rcb- This figure shows that large values of e imply SEDs from colder 
walls, farther away from the stars, which show no 10 fim silicate feature. The wall surface 
temperature varies from 136.5 i^ for the smallest Ret, and 80.7 i^ for the largest Rcb- In 
Table [2] we label these models as Series Ecc, where all the parameters are shown. 

Figure [6] shows the effects of changing the angle a, for = 0.5,cos2 = 0.5, -Rcb = l-7a and 
h = 0.28a. Results are shown for a = 0, 0.25, 0.5, 1.5, 1.75. In order to see the differences in 
terms of e, we show results for e = and e = 0.2. For symmetry arguments (with respect to 
the = axis), the configurations for a = 0.75 and 1.25 correspond to the same flux. The 
same can be said, for the a = 0.25 and 1.75 cases. Moreover, the flux for all these cases is 
similar because the stars' masses are not so different and the binary orbit is circular. The 
configurations for a = and 1 correspond to the same flux, thus, the latter is not shown. 
The minimum flux is given in the a = 0.5 model, where the disk is inclined around a line 
perpendicular to the line that connects both stars, for a = 0.5 and a = 1.5 the 10/im band 
disappears. The maximum flux occurs in the a = model (or a = 1). Changing a does not 
modify the area of the emitting region, since it only rotates this region in the plane of the 
sky. However, since the stars are not at the center of the wall, when a changes the distances 
of both stars to the visible area of the wall also change, and thus the resulting SED changes. 
We can conclude from Figure [6] that the cases are divided in two classes. The flrst one does 
not have a lOfim band and the second one has this feature. These models, where a changes, 
for cases e = and e = 2 are shown in Table [2] as Series a, ecc- 

Figure [7] shows the SED when i and h varies, for e = 0,0 = 0.5, a = 0, and Rcb = 1.7a. 
The left plot corresponds to h = 0.28a and cosi = 0.3,0.5,0.7, and the right one is for 
cosi = 0.5 and h = 0.14,0.28,0.56. The SED depends weakly on i, but, it has a strong 
dependence on h, namely, the flux increases for larger h. The reason for the former, is that 
for a perfectly vertical wall, the flux does not change monotonically with i. Of course, this 
effect will change if a realistic shape of the wall is taken. However, at this moment there are no 
studies (as far as we know) which define this shape (see Section El) ; thus, we assume a vertica l 



wall. Also, we are ignoring possible oscillation modes in the wall (JHayasaki fc Okazakill2009l ) 
due to the tidal potential of the binary. If the disk is viewed pole-on there is no visible 
area. In the opposite case, the emission from the wall of a disk seen edge-on is completely 
extinguished by the disk itself. Thus, the range in cosi relevant for modeling is small and 
the variation can be neglected. In the present model, the height of the wall is taken as an 
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arbitrary free parameter that can be adjusted to compensate for any variation of i. Models 
showing changes in i and h are presented in Table |2] as Series i,h. 

A complete analysis of all the possibilities that describe the dust composition of the wall 
is computationally demanding, so we decided to explore the changes on the SED for typical 
cases. All these models can be used as a starting point to tackle future mid-IR spectra 
modeling with observations of new binary targets. 

Figure M shows the effect on the SED of changing the type of silicates, for e = 0,0 = 
0.5, a = O,cosz = 0.5,Rcb = 1.7a and h = 0.28a. We show results for the fiducial model, and 
for the case where the silicates component are changed from pyroxenes {Mqqs Fe^^ SiO^) 
to olivines {Mg FeSiOi), model POl in Table [2] and the only member of Series pyr,oliv. 
The flux from the wall is larger when the dust is made by olivines (model POl) than when 
it is made by pyroxenes (fiducial model). This can be explained, because the former model 
has a higher temperature at the same distance from the central star. Another difference is 
that the central wavelength of the lO/im feature due to olivines is slightly shifted to longer 
wavelengths than the feature for the pyroxenes. 

Figure [9] shows the effect on the SED of changing the type of carbon, for e = 0,0 = 
0.5, a = O,cosi = 0.5,-R ^ft = 1.7a and h = 0.28a. We show results for the case where we 



take amorphous carbon (JMathis &: Whiffenlll989l ) and graphite (jWeingartner &: Drainell200ll ) 



grains. First, we compare two cases: a model with no graphite nor amorphous carbon (i.e., 
only pyroxene grains in the wall, model PAl) and a model like the fiducial one, but where 
graphite has been replaced by amorphous carbon (model PA2). We can see that model PAl 
has a lower flux than model PA2, due to a lower T^jaii, as noted in Table [2l The presence of 
amorphous carbon in model PA2 results in a larger temperature. For the same abundances, 
the absorption coefficients for graphite and amorphous carbon are similar between 0.3 and 
3 microns, with the graphite opacity slightly larger. Outside this interval, the amorphous 
carbon opacity is higher. Thus, changes in the SED (see Figure Mj are expected due to the 
differences in the opacity for both components. 

Also in Figure [HI models with variations on the abundance of graphite and amorphous 
carbon are also analyzed, for e = 0,0 = 0.5, a = O,cosi = 0.5, Rcb = 1.7a and h = 0.28a. 
Results are shown for the cases where there is no amorphous carbon and the graphite abun- 
dance is halved and doubled (models PA3 and PA4 in table [2]). In another case, the graphite 
is replaced by amorphous carbon; model PA5 corresponds to half the adopted standard value 
for amorphous carbon abundance and PA6 to the doubled value. For graphite and amor- 
phous carbon, the flux increases as ( increases. Model PAl has the lowest flux. Physically 
this implies that silicates alone are heated at a much lower temperatures than the fiducial 
model, for the adopted value of Rcb- All the PA models are assembled in Table [2] as Series 
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pyr,aniorph. 

Figure [TD] shows the effect on the SED of changing the water ice abundance, for e = 0,0 
0.5, a = O,cosi = 0.5, Ret = 1.7a and h = 0.28a. The results are for four mod els, the water 



i ce to gas mass ratios are dee = 0.0056 (model PI4, abundance proposed by iPoUack et al. 



(11994) for accretion disks), 50% (model PIS), 10% (model PI2) and 1% (model PIl) of the 



Pollack et al. value. These models are labeled in Table [2 ] as Se ries pyr,ice. We also use 



optical properties for crystalline water ice given by IWarreru (119841 ). We can see that for the 
two higher values of ice abundances, the 10/im feature disappears. The 10% case distorts this 
band significantly. Thus, for only 1% of this abundance, the SED is similar to our fiducial 
model. As we see in § [5] and § [6], this constrains the maximum allowed abundance of ice in 
the case of CoKu Tau/4. The wall temperature decreases as the ice abundance increases, 
thus, the presence of the lO/im band critically depends on T, 



wall- 



Finally, Figure [TT] shows the effect on the SED at different times, for a = O,cosi = 0.5, 
and h = 0.28a. The results are shown for the evolution of the fiducial model in half an 
orbit, for times t = x T/2 = (0,0.25,0.5,0.75, 1.0) x T/2, where T is the orbital period. 
The maximum variation in the spectra occurs between t = and t = T/2. In order to 
complete the orbit, the sequence of plots must be followed in inverse order. As expected, a 
circular binary system inside a circular CB disk shows small changes, within 10 %. Also in 
Figure [m the time variation for an e = 0.2 system is shown. The results presented are for 
t = (0, 0.2, 0.45, 0.7, 0.9) x T/2. When the stars are along a line perpendicular to the major 
axis, the wall emits its largest observable flux. The spread in fluxes is great given that in an 
eccentric system, the distances between the stars and the wall differ widely along the orbit. 
All these cases are placed under the name Series t,ecc in Table [2l 



4.2. Effects of dust in the hole 

Figure [12] shows the differences in the shape of the SED produced by the material 
within the hole. We present different values for the exponent of the mass surface density, 
p = 0, 2, 10, in the case of e = 0. In this example, the optically thin material is assumed to 
be between R = Rmin = l-4a and Rcb = 1.7a and is composed of 1 x IO^^Mq of pyroxene or 
olivine. 

The emission of material highly concentrated in the wall {p = 10) is smaller than for 
the flatter density distributions {p = and p = 2) for each type of silicates. Also, there 
is considerably more emission for the hole filled with olivines than for the hole filled by 
pyroxenes. The other noteworthy difference is that in the case of olivines, a lO/xm band is 
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more conspicuous in the SED than for pyroxenes. The reason for all this is that, for the same 
system configuration, the olivines temperature is larger than the pyroxenes' temperature. 
In addition, the band for olivine peaks at a longer wavelength ( 0.5/im greater) than the 
pyroxene band, making it look more like a peak since this wavelength is in the Wien's law 
side of the planck function. 



Figure [13] shows the effect on the SED of changing Rmin, the inner boundary of material 
in the hole. We show results for 3 cases: Rmin = l-Oa,-Rmm = l-2a and Rmin = 1.4a. As 
expected, the flux is larger for the smaller value for Rmin, for which the material lies in a 
larger region of the hole. 



Figure [M] shows the effect on the SED for different abundances of silicates. The re- 
sults ar e for the case where graphite is included with the fiducial model abundance (Cgrap = 



0.0025, 



Draine fc Led (Il984j )). The three studied cases are the standard value taken from 



Draine fc Led ( 19841 ) ((su = 0.0034) as well as half and twice this value. Note that dimin- 
ishing (sii does not imply that the flux is smaller at all wavelengths because the opacity de- 
creases, but the dust temperature increases. Changing (su and keeping (grap constant means 
that we are changing the total dust to gas mass ratio and also the relative abundances of 
silicate and graphite. In this case both the fluxes and the shape of the SED change. On the 
other hand, since the emission emerges from an optically thin region, changing the dust to 
gas mass ratio while keeping the relative abundances of each ingredient constant, will make 
the shape of the SED remain the same. In this case, the fluxes will scale proportionally with 
the dust mass inside the hole. Thus, the inferred dust mass is sensitive to the dust properties 
and composition. 



5. Application to the SED of CoKu Tau/4 

The fiducial model was made using the known properties of CoKu Tau/4, and the 
geometri cal restrictions given by the dynamical interaction betw een the binary system and 
the disk ( jPichardo et al.ll2005l l2008t lArtymowicz &: Lubowlll994l ). For the latter issue, note 
that the D0 5 models required a wall radius betw een 9 and 12AU. For instance, e = 0.8, 
Rcb ^ 3.5a dPichardoet al.ll2005h . and a ~ 8 AU Jlreland fc K rausl l2008h imphes R .h ~ 28 



AU, i.e., three times the required value in D05. Even for e = 0, iPichardo et al.l (120051 ) 
find Rcb ~ 1.9a, giving Rcb ~ 15 AU, still too big for the implicitly required temperatures 
to exp lain the observed SED. However, including the disk viscosity, lArtymowicz fc Lubow 
( I1994J ) find that for e = 0, Rcb ~ 1.7a, which corresponds to ~ 13.5 AU, still larger, but 
much closer to the wall radius estimated i n D05. In this section, we show that the value for 
Rcb given in lArtymowicz fc Lubowl ( 119941 ) results in a model of CoKu Tau/4, that fits the 
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observed SED, as shown in Figure [T5l 

In Figure one can see that for e > the wall is too cold, because the lOfim band 
does not match the observations. Even in some cases there is no band at lO/xm. Thus, a 
nearly circular binary orbit seems to be more consistent with the observed SED. We present 
a model, for a 4 Myrs system (model Ml in table [2]) in Figure [15] for CoKu Tau/4 (labeled 
Series CoKu Tau/4), where only the parameters of the wall are presented. Indeed, this 
corresponds to the fiducial model (model El) but here it is repeated for clarity. 

The spectral types and masses of the stars consistent with the observed spectrum at high 
frequencies are Ml and O.SMq and MO and O.QMq, respectively. The spectrum corrected 
using the Draine's law is better fitted by pyroxene dust grains. The higher wall temperature 
when olivines are used p roduces a flux l arger than ob se rved. Figure E] presei its the obser- 
vation s dereddened with IPrainel fcoosl ). iMathisI f ll99oh . iMcClurel (I2OO9I) and iMoneti et al. 
( I2OOII ) laws. The flux for t he observation s dereddened with the iMoneti et al.l (I2OOII ) law 
is larger than that using the iDraind (120031 ) law, the latter being taken here. Thus, olivine 
is a good option for modeling CoKu Tau/4, when using the Moneti law. A change in h 
can also increase the flux, but the depth of th e 10/xm feature changes (see Figure [7]). The 
spectrum dereddened with the iMcClurd (|2009[ ) law using Ay = 2.5 is indistinguishable from 
the data dereddened with [Drama's law and Ay = 3. Thus, the model presented here can 
be used for both cases. We use h = 0.28a {h = 0A5Rcb) as a free parameter adjusted to 
fit the observed SEDs, but in the future it would be good to have an estimated range of h 
consistent with hydrodynamical arguments. Also, the wall cannot have a larger abundance 
of water ice {Qce > 5.6 x 10^); otherwise, it would change the shape of the 10 /xm observed 
band, as described in § 14.11 and shown in Figure [TOl From Figure [9] one can conclude that 
amorphous carbon instead of graphite does not help to improve the fit, because the resulting 
lO/xm band does not match the observed SED. 

In the case of CoKu Tau/4, there is no evidence for circumstellar disks; thus, only a very 
small amount of mass should be able to cross the inner boundary of the CB disk. But even so, 
it could be enough material to produce a contribution to the mid-IR SED. Our wall emission 
model for CoKu Tau/4 shown in Figure [TS] seems to produce a reasonable fit to the observed 
SED. This fit can be improved if optically thin material in the hole should be included. The 
fiducial model requires more fiux at the right side of the lO/xm band. Thus, looking at the 
inner hole emission models in § I4.2[ we choose the one with p = 1, Rmm = l-4a and pyroxenes 
dust grains (see Figure [T2|) . A model with a lower Rmin or with graphite has a prominent 
band at lO/xm, which implies a contribution at both sides of the fiducial model band, which 
is not what is needed in this case. The surface density profile is characterized by an exponent 
p = 1 (see § l2.2p . The total mass in dust in this region is 0.006 lunar masses. The importance 
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of this final fit is that the optically thin material lies outside the Lagrangian point radius 
(i^^J, i.e. the particles are located in a dynamically allowed re gion. Some of these pa.rticles 



will b e able to accrete onto the stars due to viscosity effects ( lArtymowicz fc Lubowlll994j . 
19961 ) or due to a not well-determined interior boundary (Rmin)- However, even if this mass 



is accreted, and a standard dust-to-gas mass ratio is assumed, the n iass accretion rate ont o 
the stars is below the observational threshold, M < 10~^'^MQ/year (JMuzeroUe et al.l 120001 ). 



Figure dSl shows the model with and without material in the hole and the contribution 
from the optically thin material. Thus, the best-fit model to the emission of CoKu Tau/4 is 
the fiducial model, plus optically thin material distributed inside a geometrically thin ring 
close to the inner circumbinary wall. 



6. Conclusions 

We found a model that explains the observed mid-IR SED of CoKu Tau/4. The model 
has a binary system in a nearl y circular orbit with 4 Myrs stars old (using the evolutionary 
tracks from ISiess et al.l (120001 )) and a circumbinary disk truncated at Rcb = 1.7a ~ 12.6 
AU, with a = 7.43 AU. The height of the inner wall of the CB disk is around h/a ~ 0.3, 
for i = 60°. Also, we find that the inner hole might have a very small amount of dust 
(< 0.006 lunar masses) with a sp atial distribution in agr e emeii t with dynamical argum ents 
of the binary-d i sk int eraction (see lArtymowicz fc Lubowl (Il994j ) , iPichardo et al.l (120051 ) and 
Pichardo et al.l (120081 )). A detailed study of the amount and distribution of mass in the hole, 
can give us information about the viscously driven process that allows for some particles 
to fall into the hole when loosing angular momentu m. This phenomenon stron gly depends 
on h and a (viscosity coefficient from a-disk models. IShakura fc SunyaevI (119731 ) ). thus such 
a study can indirectly characterize these parameters. The model for the optically thin 
material in the hole can be used for arbitrary distributions of material in order to study the 
st reams. The struc t ure of the streams can be taken from the hydrodynamical simulations 
Giinther fc Kleyl ( 20021 ) and iGiinther et al.l (120041 ). However, in this preliminary study 



m 



we adopt an axisymmetric dust distribution which might represent some kind of average of 
these streams'models. We leave a more detailed approach for future work. 

Being a binary in a circular orbit, we predict a small variability (~ 10 %) of the SED 
of CoKu Tau/4 in half its orbital timescale (~ 10 years), related to the changing distances 
between the stars and the CB disk. Other sources of variability are not accounted for in the 
present model. This prediction can be tested if observations at other times are available for 
this or other binary systems. On the other hand, if the binary orbit is not circular we found 
that no model can explain the observed SED because the inner wall of the CB disk is too 
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far away from the stars to have the appropriate temperature. 

The study of the inner wall and the constraints on the dust sizes and abundances allows 
us to probe the dust even at the midplane, which is where planets must form. Note that 
Spitzer is sensitive to the dust located in the atmosphere of typical disks, not for matter 
in the midplane. The wall can be seen as a perpendicular cut through the disk. Thus, the 
details of the vertical structure are clearly exposed. In particular, the midplane is the place 
where material accumulates by set tling. Eventually, the in ass in this region is above the 
gravit ational instability th reshold ( jGoldreich fc Wardlll973l ). and planetesimals are able to 
form (JYoudin fc Shull2002l ). Thus, the study of this region is of prime importance. The 
advent of ALMA will allow us to actually see the wall. The current interferometers are able 
to get a resolution of 1", and at distances of typical regions of star f ormation, that means 
resolutions ~ 15AU. ALMA will have five times better resolution (JGuilloteau fc Dutrey 



20081 ) . thus the region of optically thin material inside the wall and the location of the wall 
can be resolved. The images produced can give us hints of the planetesimal formation process 
which will result in a complementary study to the disk structure analysis as the one presented 
here. Thus, studies of circumbinary and transitional disks, in particular the analysis of grain 
growth and settling in the wall, can place these disks in an evolutionary scheme. The dust 
composition can be traced since the formation process of the star+disk system, (due to the 
fact that muc h of the crystallization,iu st to name one process, occurs in the f irst stages of 



the f ormation, iDuUemond et al.l (120061 ) ). until the phase of very evolved stars (JGielen et al. 



20071 ). The connection between evolved stars and the current solar system is seen in binaries 
like AC Her and RU Gen. The infrared emission from the GB post-AGB disk s hows a strong 



resemblance to the spectrum of the solar system primitive comet Hale-Bopp (JGielen et al. 



20071 ). Thus, studies of dust emission of the inner part of the GB disk are important to 
improve our understanding of the bridge between young and old disks, from the formation 
of the disk to the formation of a planetary system. 

In general, the model is not restricted to any configuration nor any particular com- 
bination of stars, so it could be applied and tested with other observed binaries with cir- 
cumbinary disks. It is c lear that observations with high resolution, like the ones described by 
Ireland &: Kraud (120081 ) . in which both stars can be characterized, are very useful. In t he near 
future, other binary sy stems obser ved with Spitzer: H en 3-600 (lUchida et al.ll2004l ). St 34 
jHartmann et allboosh . HD 98800 JPurlan et al.ll2007l ). will be studied with this model. For 
these cases, the assumption was made that the binary can be approximated as an equivalent 
star with the luminosity of both. However, the results that come from this work, argue that 
this approximation should be taken with care. Finally, a more detaile d w all structure can be 



adopt ed, as the inner curved rim described by llsella fc Nattal (120051 ) or iTannirkulam et al. 



(120071 ). In these works, the hole was formed by dust sublimation and the shape of the wall 
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was given in the former by a density dependent sublimation of grains, and in the latter by 
the vertical differentiation of grain sizes. In the case of a circumbinary disk, it is not so clear 
what the shape of the disk inner wall would be, but one expects that dynamical effects of 
the stars, hydro dynamical instabilities, pressure, viscosity or even evaporation of the disk 
due to radiation from the stars, will have a major role in the final shape of the wall. The 
relevance of the unrealistic assumption of a vertical wall is not the aim of this work, however 
a full characterization requires a detailed analysis, which is beyond the scope of this paper, 
and will be addressed in a following paper. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract 
with NASA. We are indebted with Adam Kraus for the opportunity to exchange ideas with 
him. PD acknowledges a grant from CONACyT, Mexico. EN appreciates postdoctoral 
fellowships from UNAM and CONACyT, Mexico. CE was supported by the National Science 
Foundation under Award No. 0901947. 



A. APPENDIX 



The cooling/heating time is: 



2fimHaT^ ' 



tc/h = ,,_ ,^, , (Al) 



where S/j = J^"'' pdR. Note that we are interested in the area parallel to the wall 
surface. The flux is equal to aT^, assuming blackbody emission at temperature T. 

Now, a typical surface den sity for disks around T Tauri stars located at lOAU is T^d ^ 



lOgcm ^ ( jD'Alessio et al.lll998l ). and the wall height in the fiducial model is if = 0.28a 



3.34 X 10 cm. Thus, an estimate for the volumetric density is 

v 
p = ^ Ri 3 X IQ-^^gcm-^. (A2) 

The emission from the wall is optically thick, such that the emission comes from a small 
radial region. In order to give an estimate for tc/h, we take AR = H, thus S/j = S^^ = 
10gcm~'^. For the temperature T, we take T =< T^aii >= 136. 5i^, which corresponds to 
the fiducial model. Thus, the cooling/heating time is: 
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te/h ^ 2.86 X 106(EH/10(?cm-2)(T^,„/135.6)-^(/i/l.)-'s. (A3) 

The orbital time is simply 

'-- (G(M.. + M...))V. =^-°«x^°'- (A4) 

We find that for typical parameters, tc/h « torb, justifying the assumption of an in- 
stantaneous adjustment of the disk inner edge to changing stellar irradiation. 
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Fig. 1. — Schematic diagram of the parameters that characterizes the system configuration 
in space. The parameters shown are M1,M2,R1, R2,Rcb,a,(j),i,a (see text for definitions), ae 
is the distance between the center of the elhpse and the star in the focus, p corresponds to 
the pericenter, the stars configuration of closest approach. The locations of the stars are 
given with respect to the center of the ellipse, and the center of the CB disk is given with 
respect to the center of mass (CM) of the binary system. The line given by the angle a is 
the one around which the system is inclined at an angle i. The figure is not to scale. See 
text for more details. 
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Fig. 2. — Schematic diagram of an inclined system with arbitrary i. In this example, the 
system is inclined around the major axis line, a = 0. Here, dap corresponds to the distance 
between the stars in the plane of the sky. The dashed area shows the fraction of the wall 
directed towards the observer. The rest of the parameters are described in the text. The 
only difference between Figure [1] and this figure is that this one corresponds to an inclined 
system. The figure is not to scale. 
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Fig. 3. — Combined spectrum of the central stars for a 4 Myr old binary system (solid 
line). The observe d mid- IR SED of Co Ku Tau/4 is cor rected by four different exti nction 
laws: iMoneti et al.l (|2001l) (upper curve), iMcClurd (120091 ) (next curve), iDraind (120031 ) (next 
curve) andlMathid (Il990l ) (lower curve). The three last curves are very similar. IRAS data 



( jFurlan et al. 



20071 ) are represented by open squares, and the 2MASS data are crosses, circles 



and stars, corresponding to Mathis, Moneti and Draine extinction laws, respectively. The 
synthetic spectrum adopted in DOS is also plotted (dotted line). 



I 

B 

o 

U 



I 
O 



-27- 







J L 



J I I L 



J I I L 



10 



20 



30 



X(/i,m) 



Fig. 4. — SEDs of non-axisymmetrics models. The wall is an ellipse with a semiminor axis 
(b) along the line connecting both stars. The binary system has an age of 4 Myrs, with 
= 0.5,0; = O,cosi = 0.5, and h = 0.28a. The models have b = Rcb — 0.4a (solid); h = 
Rcb — 0.3a (dotted); b = Rcb — 0.2a (dashed); b = Rcb — 0.1a (dot-dashed), and b = Rcb — 0.0a 
(dot-long-dashed) . 
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Fig. 5. — SEDs of models with different binary eccentricities and an age of 4 Myrs, with 
= 0.5,0; = O,cosi = 0.5, and h = 0.28a. The parameters are given in Table[2]for each model; 
model El (solid); model E2 (pointed); model E3 (dashed); and model E4 (dot-dashed). 
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Fig. 6. — Variation of the synthetic SED with a for e = (left plot) and e = 0.2 (right 
plot), with (f) = 0.5,cosz = 0.5, and h = 0.28a. For the e = case, the models are: El, 
a = 0.0 (solid),AEl, a = 0.25 (small-dashed), AE2, a = 0.5 (long-dashed), AE3, a = 1.5 
(dot-dashed) and AE4, a = 1.75 (dot-long dashed). For the e = 0.2 case, the models are: 
AE5, a = 0.0 (sohd),AE6, a = 0.25 (small-dashed), AE7, a = 0.5 (long-dashed), AE8, 
a = 1.5 (dot-dashed) and AE9, a = 1.75 (dot-long dashed). See table [2] for more details 
about the models. 
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Fig. 7. — The left plot illustrates the dependence of the synthetic SED with inclination angle 
i, with e = 0,0 = 0.5, a = 0,Rcb = 1.7a. The plotted models are: IHl, cosi = 0.3 (solid). 
El, cosi = 0.5 (long-dashed), and IH2, cosi = 0.7 (dot-long-dashed). The right plot show 
models for various h with the same values for e,(f),a and Rcb. The models are: IH3, h = 0.14 
(sohd). El, h = 0.28, (dashed) and IH4, h = 0.56 (long-dashed). 
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Fig. 8. — Change of the SED for the fiducial model (El) for two silicates: pyroxene and 
olivines, with e = 0,0 = 0.5, a = O,cosi = 0.5, Rcb = l-7a and h = 0.28a. The solid line 
corresponds to pyroxenes (model El), and the dashed line corresponds to olivines (model 
POl). 
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Fig. 9. — Change of the SED for the fiducial model (El) for various carbon components, 
with e = 0,0 = 0.5,0; = O,cosi = 0.5,-Rc6 = l-7a and h = 0.28a. The solid line cor- 
responds to the fiducial model (El, (grap = 0.0025), the dashed line corresponds to the 
model with zero graphite abundance (PAl), and the long-dashed curve shows the model 
with amorphous carbon instead of graphite with the same abundance (PA2). Also, we 
show cases with zero amorphous carbon abundance for Cgrap = 0.00125 (PAS, dotted), and 
for Cgrap = 0.0050 (PA4,dot-dashed). Other cases with zero graphite abundance include 
Cam = 0.00125 (PA5,dot-long-dashed) and (am = 0.0050 (PA6, small-long-dashed). 
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Fig. 10. — Change of the SED for the fiducial model (El), including water ice, with e = 0,0 = 
0.5,0; = O,cosi = 0.5,-Rcb = l-7a and h = 0.28a. The plotted models are Qce = 0.000056 
(PIl,solid), (ice = 0.00056 (PI2,dashed),Cice = 0.0028 (PI3,long-dashed), and Oce = 0.0056 
(PI4,dot-dashed). 
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Fig. 11. — Time evolution of the synthetic SED of an e = (left plot) and an e = 0.2 (right 
plot) system. For the e = case, the plotted models are: PEl, t = O.OT/2 (solid), PE2, 
t = 0.25T/2 (small-dashed). El, t = 0.5T/2 (long-dashed), PES, t = 0.75T/2 (dot-dashed), 
and PE4, t = l.OT/2 (dot-long dashed). For the e = 0.2 case, the plotted models are: PE4, 
t = O.OT/2 (solid), PES, t = 0.2T/2 (small-dashed), PE6, t = 0.45T/2 (long-dashed), PE7, 
t = 0.7T/2 (dot-dashed), and PES, t = 0.9T/2 (dot-long dashed). T is the orbital period. 
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Fig. 12. — Emission of optically thin material in the disk hole of a 4 Myrs old binary system. 
Pyroxenes: p = (middle), p = 2 (upper), p = 10 (lower), all solid lines. Olivines: p = 
(middle), p = 2 (upper) ,p = 10 (lower), all dashed lines. For each case, the mass of dust in 
the hole is 1 x lO^^Mp, and it is distributed between R=1.4a and R=1.7a. 
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Fig. 13. — Emission of optically thin material in the disk hole of a 4 Myrs old binary system 
for various Rmin- Rmin = 1-Oa (solid), Rmin = l-2a (dotted) and Rmm = l-4a (dashed). For 
each case, the mass of dust in the hole is 8 x 10~^Mq, 6.2 x 10"^Mq and 4.8 x 10~^Mfr 
respectively, and it is distributed between R = Rmin and R = 1.7a. 
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Fig. 14. — Emission of optically thin material in the disk hole of a 4 Myrs old binary system 
for various silicates abundances, (su = 0.0017 (sohd), (su = 0.0034 (dotted) and (su = 0.0068 
(dashed). For each case, the mass of dust in the hole is distributed between R=1.4a and 
R=1.7a. 
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Fig. 15. — SED of model Ml, with e = 0, age= 4 Myrs (see table [2]) com pared to the S pitzer 
IRS SED of CoKu Tau/4 (points), corrected with the reddening law of IDraind ( 120031 ). We 
show the total SED (solid line), the contribution from the wall (dashed line) and from 
optically thin material in the hole (dotted line). The right plot shows the details of the 
10/im band. 



Table 1. Stellar Parameters 



spectral type M,{Mq) T,{K) R,{Rq) L,{Lq) 

Ml 0.5 3894 1.178 0.2886 

MO 0.6 4027 1.261 0.3782 
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Table 2. Models discussed in this paper'' 



Scries 


Figure 


Model 


a{AU) 
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Hn) 


a(7r) 


cosi 


Rcb{a) 


h{a) 


Dust'' 


Twall 


Ecc 
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El 


7.4S 
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0.28 




1S6.5 
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E2 
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0.2 


0.5 





0.5 


2.S5 


0.28 




110.7 
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E3 
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0.5 





0.5 
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0.28 
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E4 


18.575 
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0.5 
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0.28 
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a, ecc 
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AEl 


9.40 
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0.25 
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0.28 




125.0 
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AE2 


14.87 
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0.28 
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10S.8 
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0.28 
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AE7 
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i, h 
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pyr,oliv 
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1.7 


0.28 
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pyr,amorph 
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9 


PA6 
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1.7 


0.28 


7 


17S.8 


pyr.ice 


10 


PIl 


7.4S 





0.5 





0.5 


1.7 


0.28 


8 


136.2 




10 


PI2 


7.4S 
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0.5 
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0.28 
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1S4.0 
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7.4S 
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0.28 
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0.28 
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1S5.8 
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PE2 
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PES 
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PE4 


7.4S 





1 





0.5 


1.7 


0.28 
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Table 2 — Continued 



Series 


Figure 


Model 


a{AU) 


e 


^(tt) 


a{7T) 


cosi 


Rcb{a) 


h{a) 


Dust'' 


Twall 


CoKu Tau/4 


11 

15 


PES 

Ml 


9.29 
7.43 


0.2 



1.5 
0.5 






0.5 
0.5 


2.9 
1.7 


0.28 
0.28 


1 
1 


108.7 
136.5 



"The column titles are given as in the text. 

''The numbers in the column "Dust" refer to different adopted abundances. See Table [3] for the dust abundances. 
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Table 3. Dust abundances 



Dust model (, 



grap 



■,sil 
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■,amorp 



wall siP 
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0.0 


oliv 
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0.0 
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0.0034 


7.68E-4 


5.6E-5 


0.0 


pyr 


9 


0.0025 


0.0034 


7.68E-4 


5.6E-4 


0.0 


pyr 


10 


0.0025 


0.0034 


7.68E-4 


2.8E-3 


0.0 


pyr 


11 


0.0025 


0.0034 


7.68E-4 


0.0056 


0.0 


pyr 



to 



'Wall sil" refers to the type of silicates used in the wall. 



